Local density of states around a magnetic impurity 
in high- Tc superconductors based on the t-J model 
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The local density of states (LDOS) around a magnetic 
impurity in high- Tc superconductors is studied using the two- 
dimensional t-J model with a realistic band structure. The or- 
der parameters are determined in a self-consistent way within 
the Gutzwiller approximation and the Bogoliubov-de Gennes 
theory. In sharp contrast with the nonmagnetic impurity 
case, the LDOS near the magnetic impurity shows two res- 
onance peaks reflecting the presence of spin-dependent reso- 
nance states. It is also shown that these resonance states are 
approximately localized around the impurity. The present 
results have an large implication on the scanning tunneling 
spectroscopy observation of Bi2Sr2Ca(Cui-a;Ni[Zn]a;)208+i. 

Due to the recent development in technology, scan- 
ning tunneling spectroscopy (STS) data on high- Tc su- 
perconductors including impurities now serve to probe 
their quasiparticle properti«s as well as the nature of 
the superconducting phaseaa. Among various impuri- 
ties in high-Tc superconductors, Zn and Ni impurities are 
fundamental perturbations for the ground states because 
they are believed to be substituted for Cu in the Cu02 
plane and then strongly disrupt the surrounding elec- 
tronic structure, especially the spin configuration. Thus 
it is clear that information of quasiparticle states around 
Zn or Ni impurities is useful to reveal the mechanism or 
low-tempcrature transport properties of high-Tc super- 
conductors. 

Quite recently. Pan et al. succeeded in the 
STS observation of quasiparticle resonance states in 
Bi2Sr2CaCCui_a;Zn2;)208+5 with high spatial and energy 
resolutionEI. In the vicinity of the Zn impurities, they 
found an intense quasiparticle resonance peak at near 
zero-bias in the STS spectrum. The corresponding res- 
onance states were observed to be highly localized with 
four-fold symmetry around the Zn impurities. Theoreti- 
cally, the existence of impurity-induced bound states in 
d2.2_j,2-wave superconductors was first predicted hy Bal- 
atsky et al. using the T-matrix approximationOcl. Fur- 
ther, in our previous studies based on the t-J model with 
a nonmagnetic impurityQ-S, it was shown that the reso- 
nance states induced by the impurity actually give rise to 
a resonance peak in the local density of states (LDOS) 
which can be identified as that observed by Pan et al. 
It was also shown that the resonance states are approxi- 
mately localized around the impurity. The obtained spa- 
tial dependence of the resonance states can reproduce 



the STS result if we take into account the overlapping of 
the wavefunctions between STM tip (s-wave symmetry) 
and qwasiparticles (c?a;2_j,2-wave symmetry) in the Cu02 
planm. 

Here, we must consider the next question of the quasi- 
particle states around Ni impurities. The significant dif- 
ference between Zn and Ni impurities is that, if they 
maintain the nominal Cu^+ charge, the Zn^+ impurity 
would have a (3c?) S = configuration, while Ni^+ 
a (3c?)®, S = 1 configuration. In this framework, the 
question is how the magnetic moment of the Ni impu- 
rity affects impurity-induced quasiparticle states. Al- 
though there are some theoretical works on quasiparti- 
cle states arourid-, j^-jnagnetic impurity in dx2_y2-wave 
superconductorsliH~t3, no LDOS for direct comparison 
with STS data on Ni impurities has been obtained the- 
oretically. On the other hand, STS observation of Ni- 
doped BSCCO is only a matter in time. Therefore it is 
urgent to study the LDOS around a. Ni impurity in sim- 
ilar detail to the Zn impurity caseQ. In this letter, we 
study the quasiparticle states around a magnetic impu- 
rity in the same way as in the nonmagnetic (Zn) case. 

A Zn impurity in a Cu02 plane is considered to be a 
unitary scatterer with 5 = 0. Thus, it has been, ijeaspp^ 
ably modeled by|-aj^oint-like repulsive potentialjl3'Ej~t3 
or a vacant sitet^O in the previous theoretical works. 
On the other hand, behavior of a Ni impurity in a Cu02 
plane can be rather complicated. To take into account 
the (3d)*, 5 = 1 configuration of the Ni impurity and 
da;2_j^2-wave pairing state with a short coherence length 
in high-Tc superconductop, we start here with the model 
studied by Poilblanc et a£d. This model is the t-J model 
with a magnetic impurity, which is allowed to couple to 
its nearest-neighbor spins via an exchange coupling Jq 
only (i.e., the electron transfers onto the impurity site is 
excluded.). In order to reproduce the realistic band struc- 
ture, the next-nearest neighbor hopping term is added to 
this model. Thus the Hamiltonian of the present model 
is written as 

'H = -t Pg[4.Cj. + \i.c.)Pg + J ^i-^i 

-t' Y ^G(cLcj.+h.C.)PG-/i$]4c,, 
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in the standard notation where and mean 

the summation over nearest-neighbor and next-nearest- 
neighbor pairs, So represents impurity spin operator 
with spin-1, and + t represents the nearest- neighbor 
sites of the impurity site. The Gutzwiller's projection 
operator Pq is defined as Pq = ni(l — ni-^riii). Note 
that the four bonds connected to the impurity are ex- 
cluded in the first three terms. The last term of Eq. (|l|) 
corresponds to the coupling of the impurity to the neigh- 
boring spins of the system. When the impurity coupling 
Jo vanishes, Eq. reduces to the vacancy model which 
simulates a Zn impurity in a Cu02 plane. Throughout 
this letter, we take J/t — 0.2, t' /t — —0.4 and the hole 
doping rate S = 0.15. This set of model ]jarameters cor- 
responds to BSCCO near optimal dopingEl. 

Although Eq. (0) is a simplified model for Cu02 plane 
with a Ni impurity, it is difficult to treat even in a mean- 
field theory. Therefore, we make a further simplification 
to this model, i.e., the impurity spin So is treated as a 
fixed Ising-like spin: {Sq). Thus, the last term of Eq. (|l|) 
reduces to 



(2) 



Now the effect of the magnetism of the impurity is rep- 
resented by the effective magnetic field hcs — Jo{Sq) on 
the nearest-neighbor sites of the impurity. On the basis 
of the above simplification, we can proceed to a mean,- 
field calculation using the Gutzwiller approximatior£3, 
and then we obtain the Bogoliubov-de Gennes (BdG) 
equation and a set of self-consistent equations similar to 
that in our previous worksQ: 
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with 



V 

F*i = - E J(l - ^.,o)(l - <5,- o)A.,'5,,,+., (4) 

r 

where i + t and i + i> represent the nearest- and the 
next-nearest neighbor sites of the site i, a — ±1, and 
gt, Qs are the renormalization factors in the Gutzwiller 
approximation given by 
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The self-consistent equations are 

Q 



In the following, we assume ^iji — ^iji = ^ij and that A^j 
is a singlet pairing, i.e., A^j = Aji. Since we consider a 
well-isolated impurity, /i is fixed to the bulk value /io 
determined without impurities. 

In the present calculation, we regard the 25 x 25 square 
lattice as a unit cell of which the impurity is located at 
the center. We assume a translational symmetry of Ajj 
with respect to this unit cell. We have confirmed that this 
choice of the size of the unit cell does not have boundary 
effects and that the results for an isolated impurity can 
be simulated. Then we make use of the Fourier transform 
of the BdG equation, for which we take the number of 
the unit cells Nc = 20 x 20. We solve numerically the 
BdG equation and carry out an iteration until the self- 
consistent equations for A^ and are satisfied. Using 
(wf (fc), (fc)) and E"-{k), which are the eigenvectors 
and eigenvalues of the Fourier transformed BdG equa- 
tions, we calculate the LDOS defined by 



+ \vf{kr 5{E'^{k) + E)\, {7) 



where i represents a site, a is the index of the eigenstates, 
k is the Bloch wave number to the impurity unit cells. 

First, we compare the LDOS around the magnetic im- 
purity with one around the nonmagnetic impurity. Fig- 
ure 1 shows that the LDOS obtained on the nearest- 
neighbor site for (a) /loff = (nonmagnetic) and (b) 
^cff = 0.16t. Note that the dashed lines in Fig. 1 (a) 
and (b) represent the LDOS obtained on the site lo- 
cated at the corner of the unit cell, which reproduces the 
bulk d-wave density of states with V-shaped gap struc- 
ture. Here, the superconducting gap edges recover to 
their bulk value E = ±0.1 7t for the present parameter 
choice. Thus, we can confirm that the impurity is well- 
isolated in our calculations. In both nonmagnetic and 
magnetic cases, we find peaks near the Fermi energy, re- 
fiecting resonances caused by impurity scattering. The 
energy levels of the peaks are distinctly different in two 
cases. In the nonmagnetic case {hcs — 0), the resonance 
peak is found at slightly below zero-energy in the LDOS 
although the impurity scattering is in the unitary limit. 
When t' = 0, on the contrary, tha_cesonance peak was 
found at slightly above zero-energyulij. This difference is 
due to the change of the band structure; here the Van 
Hove singularity exists at E = — 0.085i, i.e., below the 
Fermi energy. The ratio of the resonance energy to the 
energy gap obtained here (~ 5%) shows a good agree- 
ment with the experimental result (~ 3%) in Zn-doped 
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BSCCOB. In the magnetic case, on the other hand, we 
can see that the resonance peak sphts into two. Since 
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FIG. 1. The local density of states around the (a) non- 
magnetic and (b) magnetic impurity (/icff = 0.16t). Solid lines 
are LDOS obtained on the nearest-neighbor site of the impu- 
rity. Dashed lines represent the LDOS obtained on the site 
on the corner of the unit cell. 

the obtained superconducting order parameters are real, 
the peak splitting is not due to the superconducting state 
with broken time-reversal symmetry ( e.g. the d-\-is state 
) but due to the effective magnetic field hcs. These two 
peaks correspond to the up and the down spin compo- 
nents of the quasiparticlcs. The peaks are shifted from 
the resonance energy for the nonmagnetic case because 
of the energy gain (loss) during the scattering procesiil. 

It is interesting to investigate the spatial extension of 
the impurity-induced resonance states corresponding to 
the resonance peaks in Fig. 1 (a) and (b). In Fig. 2, the 
DOS with the resonance energies are plotted as a function 
of positions around (a) the nonmagnetic and (b),(c) the 
magnetic impurity. We can see that these three states are 
approximately localized around the impurity (located on 
the center). In contrast to the system near half-fiUingEI, 
the spatial oscillating behavior of these resonance states 
is not visible because its period given by the inverse of 
the Fermi momentum is incommensurate with 2a, where 
a is the lattice spacing. 

Next we study the relation between the amplitude of 
the peak splitting and the effective magnetic field /icff. 
Figure 3(a) shows the LDOS on the nearest-neighbor site 
of the magnetic impurity for various values of hes- As h^s 
is increased, the amplitude of the peak splitting becomes 




FIG. 2. Spatial variations of the local density of states 
at the resonance energy plotted over 17 x 17 sites around the 
(a) nonmagnetic and (b),(c) magnetic impurity (/loff = 0.16t). 
The resonance energies are (a) E — — O.Olt, (b) E — 0.022t 
and (c) E = -0.032t. 

larger and the two peaks become broader and smaller. It 
is thus confirmed that the splitting is due to the effec- 
tive magnetic field. The /loff-dependence of the splitting 
amplitude D is explicitly plotted in Fig. 3(b). We ob- 
serve that D is approximately proportional to /icff when 
/loff < 0.24i. 

Here, let us discuss whether the peak splitting ob- 
tained here can be actually observed in STS experiment 
or not. In order to observe the peak splitting, its ampli- 
tude should be larger than the energy resolution in STS 
spectra. If we assume that the present model qualita- 
tively describes the effects of a Ni spin in a Cu02 plane, 
the amplitude of h^s — Jo(Sq) plays a crucial role in the 
peak splitting. However in general, it is difficult to es- 
timate {Sq) and Jo. Recently, in the spin gap state of 
the t-J model, it was shown that the Ni spin (S=l) is 
partially screened by the Cu rnmnents, resulting in an 
effective impurity spin S = 1/2113. In that paper, Jo is 
also estimated as Jo = J/2. Using these values as a ref- 
erence, we obtain hcs ~ 0.25 J = 0.05i and the expected 
splitting amplitude is D ^ 0.015i from Fig. 3(b). If we 
use J = 0.13 eV as a plausible value for Cu02 plane, 
this estimation gives D ~ 9.8 meV. Of course, we should 
note that the mean-field calculation of the t-J model may 
not give a reliable numbers. However, even taking this 
point into consideration, we believe that the value of D 
is still in order of meV. In the recent STS experiments, 
sub-meV resolution in energy is already accessiblea. Thus 
the peak splitting would be observed in STS experiments 
with high resolution. In order to make the splitting am- 
plitude D ( ie. {Sq)) larger, the observation in an ap- 
plied magnetic field would be effective. At this time, it is 
a natural question whether the splitting observed under 



3 



heff 


= o.ost 


heft 


= 0.16t 


heff 


= 0.24t 


heff 


= 0.52/ 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 
E/t 




FIG. 3. (a) The local density of states obtained on the 
nearest-neighbor site of the impurity for various values of 
hcH = Jo{Sq). The dashed line represents the LDOS ob- 
tained on the corner of the unit cell, (b) The /leff dependence 
of the amplitude of the peak splitting D obtained on the near- 
est-neighbor site of the impurity. 



the applied magnetic field is due to the magnetic field 
itself. We speculate that the amplitude of Zeeman split- 
ting due to the external magnetic field of a few Tesla is 
apparently smaller than that of the peak splitting ob- 
tained here. Thus, the effect of magnetic impurity will 
be well identified by checking the difference in Ni-doped 
BSCCO and Zn-doped BSCCO under the same magnetic 
field. 

In our calculation, so far, an antiferromagnetic order 
parameter and a d^y-yfave superconducting order param- 
eter which could be locally induced around the impurity 
have not been taken into account. In particular, it is an 
interesting problem to reveal whether or not a dxy-wave 
order parameter with broken time-msjErsal symmetry is 
induced around a magnetic impurity£3tj. However in the 
standard t-J model, a dxy-wave order parameter is not 
favored, so that some extension of the model is needed. 

In summary, we have investigated the LDOS around 
the magnetic impurity in the (i^2_j,2-wave superconduct- 
ing state with short coherence length based on the t-J 
model. Within the Bogoliubov-de Gennes theory and 
the Gutzwiller approximation, we predict the splitting of 
the resonance peak in LDOS near the impurity which is 
ready to be checked experimentally. 
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